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Abstract. The validity of an approximation formula for European option prices 
under a general stochastic volatility model is proved in the light of the Edgeworth 
expansion for ergodic diffusions. The asymptotic expansion is around the Black- 
Scholes price and is uniform in bounded payoff functions. The result provides a 
validation of an existing singular perturbation expansion formula for the fast mean 
reverting stochastic volatility model. 



1. Introduction 

In the last decade, many results on asymptotic expansions of option prices for 
stochastic volatility models appeared in the literature. Such an expansion formula 
gives an approximation to theoretical price of option and sheds light to the shape 
of theoretical implied volatility surface. See e.g., Gatheral (T2| for a practical guide. 
The primary objective of this article is not to introduce a new expansion formula 
but to prove the validity of an existing one which was introduced by Fouque et 
al. f7\. We suppose that the log price process Z satisfies the stochastic differential 
equation 

piXt)dWj + ^1 - piXtfdWf 
dXf = b{Xt)dt + c{Xt)dWl 
under a risk-neutral probability measure, where (W^, W^) is a 2-dimensional stan- 
dard Brownian motion, r = [rt] stands for interest rate and is assumed to be 
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D - exp i - I Tsds 



deterministic, and b, c, cp, p are Borel functions with |p| < 1. Under mild conditions 
on the ergodicity of X, we validate an approximation 

(2) DE[/(Zr)] ^ DE[(1 + p{N)) f{Zo - log(D) - E/2 + VZN)] 

for every bounded Borel function /, where N ~ N{0,1), E = n[<p^]T with the 
ergodic distribution H of X and 

Jo 

(3) p(z) = a|l-z2 + ^(z3-3z) 

I Ve 

In particular, we have a simple formula 

DE[(X - exp(Zr))+] ~ P^siK E) - adiiK L)DK(p{d2{K, E)) 

for put option price with strike K, where Pbs(^^/ ^) is the Black-Scholes price of the 
put option 

Pbs(X, E) = DK0i-d2{K E)) - exp(Zo)0(-d2(-K, E) - Ve), 

log(K)-Zo+log(D) Ve 
d2(i^,E) = — — . 

Ve 2 

Notice that if a = then the right hand side of ^ coincides with the Black-Scholes 
price for the European payoff function / o log with volatility n[(p2]^^2. The term 
with p is small if c is large, so that in such a case it should be regarded as a cor- 
rection term to the Black-Scholes approximation. The right hand side of l|2) is an 
alternative representation of the so-called fast mean reverting or singular perturba- 
tion expansion formula and its validity has been discussed by Fouque et al. ||8ll9l , 
Conlon and Sullivan [4] and Khasminskii and Yin L14J under restrictive conditions 
on the payoff function / or on the coefficients of the stochastic differential equation 
(H). Recently, Fukasawa flTl gave a general framework based on Yoshida's the- 
ory of martingale expansion to prove the validity of such an asymptotic expansion 
around the Black-Scholes price for a general stochastic volatility model with jumps, 
which in particular incorporates the fast mean reverting case with This paper, 
on the other hand, concentrates on the particular standard model to improve the 
preceding results mainly in the following points: 

(1) conditions on the integrability of (Z) are weakened, 

(2) precise order estimate of approximation error is given. 
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The framework of Fukasawa [TT) is too general to give such a precise estimate of 
order of error. A PDE approach taken by Fouque et al. ||8)|9) gave order estimates 
which depend on the regularity of the payoff /. The order given in this article is 
more precise and does not depend on the regularity of /. We require no condition 
on the smoothness of / and a weaker condition on the coefficients (p, p, b and c. We 
exploit Edgeworth expansion for ergodic diffusions developed by Fukasawa [lOJ. 

The Edgeworth expansion is a refinement of the central limit theorem and has 
played an important role in statistics. There are three approaches to validate the 
Edgeworth expansion for ergodic continuous-time processes. Global(martingale) 
and local(mixing) approaches which were developed by Yoshida fT9l and |20| re- 
spectively are widely applicable to general continuous-time processes. The third 
approach, which is called regenerative approach and was developed by Fuka- 
sawa IfTOl extending Malinovskii 1151 , is applicable only to strong Markov processes 
but requires weaker conditions of ergodicity and integrability. The martingale ap- 
proach was applied to the validation problem of perturbation expansions by Fuka- 
sawa fTTI as noted above. The present article exploits the regenerative approach 
which enables us to treat such an ergodic diffusion X that is not geometrically 
mixing. An extension to this direction is important because empirical studies such 
as Andersen et al. Ill showed that the volatility process appears "very slowly mean 
reverting", that is, the autocorrelation function decays slowly. Our model ^ with 
a condition of ergodicity given below is a natural extension of the fast mean re- 
verting model of Fouque et al. Q (§1 but does not necessarily imply a fast decay of 
the autocorrelation function. It admits a pol5momial decay of a-mixing coefficient. 

It should be noted that our approach in this article utilizes the fact that X 
is one-dimensional in ([ij. See Fukasawa ITTl for multi-dimensional fast mean 
reverting stochastic volatility model with jumps. In Section 2, we review the fast 
mean reverting expansion technique. The main result is stated in Section 3 with 
examples. An introduction to the Edgeworth expansion theory is given in Section 4 
and then, the proof of the main result is presented in Section 5. The proof of an 
important lemma is deferred to Section 6. 
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2. Fast mean reverting stochastic volatility 

2.1. PDE approach. Here we review an asjmiptotic method introduced by Fouque 
et al. IZil/ where a family of the stochastic volatility models 
dS'l = rS;'df + <p(X;')S;'dWf , 

It t J '7 

is considered, where W = (Wf) and W = (W^'') are standard Brownian motions 
with correlation (W, WP)t = pt, p e [-1, 1]. This is a special case of with p(x) s p, 
t'(x) = (m - x)/if - V V2A(x)/7], c(x) = v V2/j], where m, v are constants and A is 
a Borel function associated with the market price of volatility risk. For a given 
payoff fimction / and maturity T, the European option price at time t < T defined 
as 

(5) PHt, s, v) = e-''^^-^^E[f{S'i.)\S'l = s, X'^ = v] 

satisfies 



where 



\£o + -£i + XzlP'' = 0, PHT,s,v) = f{s) 
Xo = v— + 

£i = yl2pvs(p{v)^—- - yl2vA{v)—, 
dsdv dv 

d 1 , s2 2 '^^ / ^ IN 

Notice that Xo is the infinitesimal generator of the OU process 

(6) dX° = {m- X°)di + v V2dWf 

and X.2 is the Black-Scholes operator with volatility level |^(c)|. By formally ex- 
panding P'' in terms of /] and equating the same order terms of /] in the PDE, one 
obtains 

(7) P'' = Po + ?]Pi + higher order terms of 7] 

for the Black-Scholes price Pq with constant volatility no[^^]^^^, where Hq is the 
ergodic distribution of the OU process X°, and 

(8) Po + nPr =Po-{T-t) (V2s2^ + Vss^^ 
with constants V2 and V3 which are of 0(7]). 
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As a practical application, Fouque et al. [ZJ proposed its use in calibration 
problem. They derived an expansion of the Black-Scholes implied volatility 0bs of 
the form 

from where K. is the strike price, S is the spot price, T - Hs the time to the 
maturity, a and h are constants connecting to V2 and as 

(10) V2 = a{{a-b)-a{r+^a'')), V^ = -aa\ = Uo[(p'\. 

The calibration methodology consists of (i) estimation of from historical stock 
returns, (ii) estimation of a and h by fitting ||9| to the implied volatility surface, and 
(iii) pricing or hedging by using estimated a, a and h via l|8]l and i fTOl l. This approach 
captures the volatility skew as well as the term structure. It enables us to calibrate 
fast and stably due to parsimony of parameters; we have no more need to specify 
all the parameters in the underlying stochastic volatility model. The first step (i) 
can be eliminated because the number of essential parameters is 2 in light of lO; 
by using n,j[(p^]^^^ instead of no[^^]^^^ for a, where n,j is the ergodic distribution 
of X'', we can see that the right hand side of lHJ coincides with that of ^ with 
Va = -an,,[^2] V2 = 2V3. 

It should be explained what is the intuition of rj —> 0. To fix ideas, let A = for 
brevity. Then Xf := X'j^^t satisfies 

dXf = {m- Xt)dt + V V2dWf, 
where Wt - rf^ W,j2t is a standard Brownian motion, and it holds 

dS|' = rS'jdt + (p(X,/,,2)S|'dWf . 
Hence r] stands for the volatility time scale. Note that 

<log(S'')),= /(X,/,,2)2ds ~ rf <p(X°)Ms ^ no[<p2]i 
Jo Jo 

by the law of large numbers for ergodic diffusions, where X" is a solution of 
This convergence implies that the log price log(S|') is asymptotically normally 
distributed with mean rt - Tlo[(p^]t/2 and variance no[(p^]f by martingale central 
limit theorem. The limit is nothing but the Black-Scholes model with volatility 
no[(p^]^^^. The asymptotic expansion formula around the Black-Scholes price can 
be therefore regarded as a refinement of a normal approximation based on the 
central limit theorem for ergodic diffusions. 
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2.2. Martingale expansion. Note that a formal calculation as in (0 does not ensure 
in general that the asymptotic expansion formula is actually valid. A rigorous 
validation is not easy if the payoff / or a coefficient of the stochastic differential 
equation is not smooth. See e.g. Fouque et al. [8J. A general result on the validity 
is given by Fukasawa flT]. Here we state a simplified version of it. Consider a 
sequence of models of type ^ : 
1 



p(X«)dWi + ^1 - p(X«)2dWf 



dx;' = bn{Xl')dt + c„{Xl')dWj, 
where bn and c,„ n e N are sequences of Borel functions. 

Theorem 2.1. Suppose that for any p > 0, the LP moments of 

(11) j\{Xlfdt, !^j\{X'^f{l-p{Xlf)dt^ 

are bounded m w e N and that there exist positive sequences e„, E„ with e„ — > 0, 
Zoo := lim„^oo > such that 

(Ml (M")t-E„\ 

(12) h^ / \ " \^N{0,V) 

in law with a 2x2 variance matrix V = [Vjj] asn ^ oo^ where M" is the local martingale 
part ofZ. Then, for every Borel function f of polynomial growth, 



(13) IE[/(Z^)] = E[(l + p„(N))/(Zo - log(D) - E„/2 + V^nN)] + o{e„) 
fls n — > oo, where N ~ A/'(0, 1), D ;s defined as in (|3| and 

p„(z) = {- V^(z' - 1) + (z' - 3z)} . 

An appealing point of this theorem is that it gives a validation of not only 
the singular perturbation but also regular perturbation expansions including the 
so-called small vol-of-vol expansion. It is also noteworthy that the asymptotic 
skewness Vu appeared in the expansion formula is represented as the asymptotic 
covariance between the log price and the integrated volatility. Our interest here is 
however to deal with the singular case only. Now, suppose that b„ and (1 + c^,)/c„ 
are locally integrable and locally bounded on ]R respectively for each n e N; we 
take R as the state space of X" by a suitable scale transformation. Further, we 
assume that s„(R) = ]R for each n e N, which ensures that there exists a unique 
weak solution of l(l). See e.g., Skorokhod lllTl . Section 3.1. It is also known that the 



EDGEWORTH EXPANSION 



7 



ergodic distribution Tl„ of X" is, if exists, given by 

n„(dx)= ^ J^"" s„(x)= r exp(-2 r ^f^dwldv 
els'„{x)cl{x) Jo I Jo c«(ii')2 J 

with a normalizing constant e^: 

J s;,(a:)c2(x)' 

Theorem 2.2. Suppose that 

(1) /or any p > 0, the LP boundedness of the sequences (III] > /zoZds, 

(2) e„ ^ fls M ^ oo, 

(3) lim„^oo Tln[(p^] exists and is positive, 

(4) lim„^oo n„[(ppi/'„] and lim„^oo Tlniipl] exist, ivhere 

^n{x) = 2e„c„{x)s'„{x) I ((p(;])2 - n„[(p^])n„{di]), 

— oo 

(5) the sequences 

JX'' ^-nV-) 1 Jo 

and 

^ ./■„(X«)p(X«)<p(X[')df-n„[i/;„p(p] 
converge to m probability asn ^ oo. 

Then, the approximation ^ is valid in that ( f73l ) holds with E„ = n„[(p^]r and Pn = P 
defined as Q with b = b„,c = c„ and L = L,,. 

Proof: Let us verify ^12) with E„ = n„[(p^]r and 

Vi2 = -2 lim E-3/2j-[„[<ppi/;„]. 
Notice that by the Ito-Tanaka formula, 

<M")r - n„[(p2]r = Cicpixif - n„[(p2])df 

Jo 

= e„ ^dx-e„ i/'„(X|')dWi. 
Jx;' c„(,xj Jo 

It suffices then to prove the asymptotic normality of 

(P(x«) [p(xndwi + ^1 - p(x;')2dwf] , !/^„(x|')dwiJ. 

This follows from the martingale central limit theorem imder the fifth assumption. 

//// 
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The conditions are easily verified in such a case that both Il„ and s„ do not 
depend on n e N. The model l|4| with A = and /] - i]„, where rj„ is a positive 
sequence with i]„ 0, is an example of such an easy case. 

3. Main results 

3.1. Main theorem and remarks. Here we state the main results of this article. We 
treat ^ with Borel functions (p, p satisfying \p\ < 1, b being a locally integrable 
function on R and c being a positive Borel function such that (1 + c^)/c is locally 
bounded on R. We suppose that (p also is locally boimded on R and that there 
exists a non-empty open set U c R such that it holds on U that 

(1) (p and p are continuously differentiable, 

(2) (1 - p2)^2 > and \(p'\ > 0. 

If (p is constant, then the approximation ^ is trivially valid. Since U can be any 
open set as long as it is not empty, this condition is not restrictive in the context 
of stochastic volatility models. This rules out, however, the case p = 1. We can 
introduce alternative framework to include such a case although we do not go 
to the details in this article for the sake of brevity. We fix (p, p, U and assume 
(Zo, Xo) - (0, 0) without loss of generality. 

Define the scale function s : R ^ R and the normalized speed measure density 
71 : R ^ R as 

with 

^^^^ ''^Jv^cjc^y 

It is well-known that the stochastic differential equation for X in ([ij has a unique 
weak solution which is ergodic if e < oo and s(R) = R. The ergodic distribution 
n of X is given by n(dx) = n{x)dx. See e.g., Skorokhod llTll . Section 3.1. Notice 
that X is completely characterized by {n,s,e). In fact, we can recover b and c by 
1/c^ = e^s'n and b = -c^s" /2s'. Taking this into mind, denote by C the set of all 
triplets (tt, s, e) with n being a locally bounded probability density function on R 
such that I/tt is also locally bounded on R, s being a bijection from R to R such that 
s' exists and is a positive absolutely continuous function, and e being a positive 
finite constant. 
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For given y = (y+,y_) e [0, oo)^ and 6 e (0,1), denote by C(y, 6) the set of 
6 = {n,s,e) e C satisfying Conditions 13. 1[ l3^ below. 

Condition 3.1. It holds that 

(1 + (p{x)^)n{x)s'{y) < exp{- log(6) + y+x - (4y+ + 6){x - y)} 
for allx> y >0 and 

(1 + q){xf)n{x)s' {y) < exp{- log(6) - y-X + (4y_ + b){x - y)] 
for all X < y < 0. 

Condition 3.2. There exist x e U and a e[6,l /d] such that \x\ < l/6,[x- a,x + a] (Z U, 
71 is absolutely continuous on[x- a,x + a] and it holds 



sup 

ife[x-a,x+a] 



7^<pp] (y) 



Vs'(y)V.(y)V-J-^V^<l/6. 



Given 6 e C, we write ng, sg, eg, bg, cg, for the elements of = {n, s, e), the 
corresponding coefficients b, c of the stochastic differential equations, and the log 
price process Z defined as ((l) respectively. 

Theorem 3.3. Fix y = {y+, y_) e [0, oo)^ and 5 e (0, 1). Denote by the set of the Borel 
functions bounded by 1/6. Then, 

sup E[/(Z?)] - E[(l + pe{N))fi- log(D) - Lg/2 + V^^N)] 

/eSs,eeC(y,6) 

is finite, where N ~ N{0,1), Eg = ng[cp^]T,ng{dx) = ng{x)dxandD,p = pg aredefined 
by ^ with E = Eg, n = Hg, c = Cg. 

Remark 3.4. The point of the definition oiC{y, S) is that it is written independently 
of e. As a result, if e C(y, 6), then (ti,|, s^,, e,j) associated with the drift coefficient 
brj = bg/ if- and the diffusion coefficient c,j = cg/rj is also an element of C{y, 6) for 
any rj > 0. In fact 7i,j = ng and s,j - sg. On the other hand, e,j - i]eg, so that 
Theorem l3.3l implies, with a slight abuse of notation, 

(16) ElfiZl)] = E[(l + p,,(N))/(Zo - log(D) - E,,/2 + ^N)] + 0(rf) 

as 7] ^ 0. 

Remark 3.5. Given e C, Condition 13.21 does not hold for any 6 > only when 
considering vicious examples such as the case ( yJng/s'gCppY is not continuous at 
any point of Lf; a sufficient condition for Condition 13.21 to hold with some 6 > 
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is that ( ^ne/s'gCppy is continuous at some point of U. If Condition |3j2] holds with 
some 6 > 0, then it holds with any 5 e (0, 6] as well. 

3.2. Examples. 

Lemma 3.6. Let 6 eC. If there exist {y+, y-) e [0, oo)^ such that 

1 + (piv)^ 

(17) K± > 2y±, km sup , , < oo 

i,^±oo e' ±l''lce(57)2 

K+ = - lim sup — -— , K- = lim mf — -— , 

j,^^ ce{vy f^-oo ce(i?)'^ 

then there exists do > such that for any 6 e (0, 6n Al), Condition \3.1\ holdsfor 6 = {■n,s,e) 
with y = (y+,y_) and 6. 

Proof This is shown in a straightforward manner by lHHl. //// 

Example 3.7. Consider 

dZ, = {n - ^Vt^dt + ^|Vt{pdWj + ^1 - p2dWf )) 

dy, = ^7]-2(f( - y^dt + /]-Vtr'dwi 

for positive constants <5, /.i, i] > 0, p e (-1, 1) and v e [1/2, oo). We assume ^p > 1/2 
if V = 1/2. Then, the scale fimction of V satisfies s^((0, oo)) = ]R, so that we can 
apply Ito's formula to X = log(y) to have 

dXf = 7f^{^pe-^' -£,- e-2(i-^')^'/2)df + rT^e-^^-^^^'dWl. 

In this scale, (p{x) - exp(x/2), so that we can take any open set as LI c ]R. We fix 

£,,p,v,p arbitrarily. In the light of Remark 13.41 it suffices to verify Conditions 13.11 

and 13.21 only when /] = 1. It is trivial that Condition 13.21 holds with a sufficiently 

small 6 > 0. If V = 1/2, then lUZll also holds with 

1 

K+ = oo, K- = E.p - -, y+ =2, y_ = 0. 
If V e (1 /2, 1), then it holds with 

K+ = oo, y^ - 3 - 2v, y- = 0. 

If V - 1, it then holds with 

1 

K+ + -, K- = oo, y+ = 1, )/_ = 

provided that > 3/2. Unf ortimately, iflTll does not hold if v e (1,11/8]. If v > 11/8, 
it then holds with 

1 

K+ - 2' ^ ^ y+ = (3 - 2v)+, y- - 2v - 2. 
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Note that the case v = 1/2 corresponds to the Heston model. In this case, we have 
a more explicit expression of the asymptotic expansion formula; we have | [T6l l with 

This is due to the fact that the ergodic distribution of the CIR process is a gamma 
distribution. 

Example 3.8. Here we treat l|4j. In order to prove the validity of the singular 
expansion in the form ((T6|l for Zj = log(S!J.), it suffices to show that there exist y, 
5 and tjq > such that Conditions 13.11 and 13.21 hold for 6 - {n,s,e) e C associated 
with 

be{x) -m-x-i-jv V2A(x), Cg{x) = v V2 

for any j] e (0, rjo], in the light of Remark I3l4l Here we fix m e ]R and v e (0, oo). 
Suppose that there exists y_) e [0, oo)^ such that 

limsupe~'^-''^''(p^(x) < oo 

and that A is locally bounded on IR with 

Aoo limmf > -oo. 

U'Hcx) X 

Then we have 

be{v) 



-sgn{v)- ■ 
cg{vy 

as \v\ oo uniformly in r] e (0, 7]o] with, say, rjo = 1 A |l/(2vAo<, A 0)|. Hence, there 
exists 6 e (0,1) such that Condition 13.11 holds for any /] e (0, r]o] with y = (y+, y_) 
and 6. By, if necessary, replacing (6, )jo) with a smaller one. Condition 13.21 also is 
verified for any i] e (0, r/o] under a slight condition on (p stated in the beginning 
of this section. Consequently, by Theorem 13.31 we have i fTsl l for iSJl if |p| < 1 
in addition. The obtained estimate of error 0{Tf) is a stronger result than one 
obtained by Fouque et al. I8H9J . 

Example 3.9. Here we treat a diffusion which is not geometrically mixing. Consider 
the stochastic differential equation 

f]^ V2 / cosh(yf)^ f]cosh(Xf) 

with <5 > 1/2 and rj > 0. Putting Yf = sinh(Xf), we have 

1 cy 1 
dY, = ---^di+-dWt 
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This stochastic differential equation has a unique weak solution which is ergodic. 
A polynomial lower bound for the a mixing coefficient is given in Veretennikov [ 18 1 
which implies in particular that X = sinh~^(Y) is not geometrically mixing for any 
4. Now, let us verify Conditions 13. ll and l3.2l for with 

1 /I \ tanh(x) , , 1 1 



r/^ \2 / cosh(x)^ 77 ( 



cosh(x)2 7] cosh(x) 

for any r] > 0. In the light of Remark l3!4l it suffices to deal with the case ?] = 1- 
Since 

- lim sgn(x)— — = +6 

|.vHoo " c{xY 2 

we have i flTt if there exists > such that 

2 



\ 

sup e"f''-^'(p(x)^ < 00, - + £ > 4 + 2/i. 



Condition 13.21 also is satisfied with a sufficiently small 6 > if the condition stated 
at the beginning of this section holds. 

4. Edgeworth expansion 

4.1. Gram-Charlier expansion. Here we give a brief introduction to the Edge- 
worth expansion. It is in a sense a rearrangement of the Gram-Charlier expansion. 
Let Y be a random variable with ]E[Y] = and E[Y^] = 1. If it has a density py with 
an integrability condition 

(18) J py(z)'</'(z)"'dz < cx,, 

where (\> is the standard normal density, then we have 



i=0 

in L^((/)) with Hermite polynomials Hj defined as the coefficients of the Taylor series 

00 ■ 

(19) e*''-*'^^ = Y^Hj{x)-, {t,x)eW}. 

j=o J- 

This is an orthonormal series expansion of pz/(p ^ ^^i'P) ^^^^ implies that 

(20) ]E[/(Y)] = -E[H,(Y)] f{z)Hj(z)(p{z)dz 

for / e L^{(p)- The Edgeworth formula is obtained by rearranging this Gram- 
Charlier series. For example, if Y = m~^^^ TJJ=i with an iid sequence Xy, then the 
7-th cumulant of Y is of 0{m^~i^^). This is simply because 
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where \pY and i/'x are the characteristic functions of Y and Xj respectively. Even if Y 
is not an iid sum, kJ = 0{m^~i^^) often remains true in cases where Y converges in 
law to a normal distribution as m ^ co. Since E[Ho(y)] = l,E[Hi(y)] = E[H2(y)] = 
and for ; > 3, 

]E[H,(y)] = Y, E 

k=l I'lH — t-rt=j, rj>3 

by (HD, it follows from (|20j that 



rj! . . . f/t! k\ 



IE[/(y)] = mqj{z)<p{z)dz + oi 



j=o 

with suitable polynomials qj. Taking / = 0, we have the central limit theorem; in 
this sense, the Edgeworth expansion is a refinement of the central limit theorem. 
This asymptotic expansion can be validated under weaker conditions than l fT8)l ; 
see Bhattacharya and Rao Q and Hall ||T3l for iid cases. Here we give one of the 
validity theorems. 

Theorem 4.1. Let X" be a triangular array of d-dimensional independent random variables 
with mean 0. Assume that X" ~ Xj for all j and that 

supE[|X;'|^] < CO 

for an integer ^ > 4, 

sup \^"{u)\ < 1, sup 1 \^"{u)\'^du < CO, 

for allb > and for some rj> 1 respectively, where 

V"(m) = E[exp(m ■ X!^'}]. 

Then, there exists mo such that SJJ, = TJjLi X" has a bounded density pjj, for all 
m > Mq, n e N. Further, it holds that 

sup ot(1 + \xf)\pl{x) - ql{x)\ < CO, 

3r6K'',m>mo,)ieN 

where 

d 



q"Jx) = cp{x;0,v")--^ Y K"d,djdk(P(x;0,v") 
with the variance matrix v" ofX" and the third moments k".^, ofX". 



Proof: This result is a variant of Theorem 19.2 of Bhattacharya and Rao (2). 
Although the distribution of X^' depends on n, the assertion is proved in a similar 
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manner with the aid of Theorem 9.10 of Bhattacharya and Rao |2]|, due to our 
assumptions. For example, we have ( and use ) 

0< inf E[|M-X^p]< sup E[|m -Xi'p] < oo. 

//// 



l"l=l,neN |u|=l,«eN 



4.2. Edgeworth expansion for regenerative functionals. We have seen that the 
fast mean reverting expansion gives a correction term to the Black-Scholes price 
that corresponds to the central limit of an additive fimctional of ergodic diffusion 
in Section 2.1. In order to prove the validity of the expansion, it is therefore natural 
to apply the Edgeworth expansion theory for ergodic diffusions. Here we present 
a general result for triangular arrays of regenerative functionals, which extends 
a result for additive functionals of ergodic diffusions given by Fukasawa [10|. 
Let P" = {D,",T'", {F"},P") be a family of filtered probability spaces satisfying the 
usual assumptions and K" = (X") be an {F"}-adapted cadlag process defined on 
P". Denote by £"[■] and Var"[ ] the expectation and variance with respect to P" 
respectively. For a given sequence of increasing (F"}-stopping times {x"} with 
- and lim^^oo x" - oo, put 

We say that K" is a regenerative functional if there exists {t"} such that 

(i) CKp rp is independent to F^„ for each ; = 1, 2, . . . , 

(ii) CTC", rp, ;• = 1, 2, . . . are identically distributed. 

An additive functional of an ergodic diffusion is a regenerative functional. See 
Fukasawa flOl for the details. Let K" be a i-dimensional regenerative functional 
and put 'k'^ = (TC",,,, ip for / = 0, 1, . . . . Notice that IC'', ; > 1 is an iid sequence 
and independent of 'Kg. Assume that Yar"['KJ] exists and is of rank d' + 1 with 
1 < d' < d for all / > 1- Without loss of generality, assume that there exists a 
rf' -dimensional iid sequence G", ; > 1 such that the variance matrix of (G", /") is of 
full rank and that 

(21) ICJ^iGpRpip 

with ad - d' dimensional sequence R". Put 

ml = E"[l\l ml = E"[G'll ;< = E"[R^], 

^« = (^;;) = K, 
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and 

= {Oj, ip, G'l = G] - \]mllml, j e N. 
Due to the definition, it is not difficult to see a law of large numbers holds: 

in probability as T ^ oo. Further, a central limit theorem 

VT(j<:«/r-fi")^//(o, V") 

holds with a suitable matrix V" . Our aim here is to give a refinement of this central 
limit theorem. More precisely, for a given function A" : R'' ^ ]R and a positive 
sequence T„ —> oo, we present a valid approximation of the distribution of 

^/T„(A"mjT„)-A"{^")) 

up to 0{T~^) as n —> CO. As far as considering this form, we can assume without 
loss of generality that E"[|R"|] = for all / > 1 in Put 

i^il,) = Var«[G;']/<, p" = ipl) = Cov"[Gl,q] 

and 

..n /^M ^n,,n ^n,,n ,,n \ i ..^ji 

f^Um = (T^Um - Pk^m " Pi ^m.k ' Pmh,i>/'^L' 

where (Kf , ) is the third moments of G". 

^ k.l.m' 1 



Condition 4.2. It holds that 



inf m" > 0. 



neN 



Condition 4.3. For £, = {d' + 2)V 4, it holds that 



sup E"[|-7Co"n+E"[|Kjn + E« 

neN 



Jt" 



dt 



< oo. 



Under Conditions 14.21 and 14.31 the sequences n", ij^li), {l^li„) are bounded in 
neN. 



Condition 4.4. Let W" be the characteristic function o/K"; 

W'iu) = E"[exp{!M ■ K'l)]. 

It holds 

sup \W"{u)\<l 

\u\>b,n€K 

for allb > and there exists rj >1 such that 



sup r \^"{u 

neN Jk'''+i 



)|''dM < CO. 
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Note that under C onditions I4.3l and l4.4[ it holds 

< inf E"[\a ■ K"^\^] < sup E"[\a ■ K'^f ] < co, 

M='^,n<^^ |fl|=l,KeN 

that is, the largest and smallest eigenvalues of the variance matrix of K" is bounded 
and bounded away from in n e N. 

Let B„(C) = \xe ]R''; \x - < Q for C > 0, 

flf = <9,A«(f("), a'l- = d,djA'\ii"), l<i,j<d 

for a given function A" : M!^ —> M. which is twice differentiable at the point fi" and 



kj=l 



Condition 4.5. There exists C,> such that 

(1) A" -.M!^ ^M. is four times continuously differentiable on B„{C,)for all n, 

(2) all the derivatives up to fourth order are bounded on B„(C) uniformly in n, 

(3) it holds that 

< inf v" < sup v" < oo. 



neN 



iieN 



sup T„ 



< oo. 



Denote by t the natural inclusion: M!^' 3 v {v,0,. . .,0) e'R.'^. 

Theorem 4.6. Let M be a positive constant and Sm be the set of Borel functions on M 
which are bounded by M. Under Conditions \4.2\\4.3\\4.4\ and \4.5\ it holds that 

E«[H( ^|f;,iA"mJTn) - A"i^"m - J H{z)q"{z)dz 
where q" is defined as 
(22) 

with <p{z; c") being the normal density with mean and variance v", 
qi{z;v") = -dcp{z;v"), qsiz-v") = -<9>(z;z;"). 



^3 



q"{z) = c/)(z;z;«) + T;,''' \A",qi{z;v") + -fq3{z;v") 



and 



(23) 



t(p") 



m. 



k,l,in=l 



j,k,l,m=l 



Proof: The proof is a repetition of the proof of Theorem 4.1 of Fukasawa lITOl 
with the aid of Theorem l4.1l and so is omitted. //// 



edgeworth expansion 
5. Proof of Theorem 
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Here we give the proof of Theorem 13.31 We are considering ^ with (p, p, 
U satisfying the condition stated at the beginning of Section 3. The initial value 
(Zo, Xo) = (0, 0) and the time to maturity T are fixed. Now, to obtain a contradiction, 
let us suppose that the supremum in Theorem 13.31 is infinite. Then there exists a 
sequence d„ e C{y, 6) such that 

^-2 



(24) 



E[/(Z!f)] - E[(l + p„(N))/(- log(D) -LJ2+ ^JL„N)] 



asn ^ CO, where e„ - eg,,, p„ = pg,,, E„ - Eg,,. Put bn - e^bg^^, c„ - e„cg^^ and denote 
by E" the expectation operator with respect to the law of X" determined by the 
stochastic differential equation 

dXf = b„{X^)dt + c„(X«)dWf\ X'^ ^ x e R, 

where is a standard Brownian motion. It is easy to see that the law of X = {Xf} 
in ^ is the same as that of < X" ^ [ with X" = 0. Hence, the law of Z" = Z^" is the 
same as that of 

under EJJ, where 

(25) Tn = \, A"ix, y) = h„ = Tcp" - E„, 

x;' = |£/;„(x::)ds,£<p(XD p(XDdwi + ^1 - p{x«Ydwl 

and (W^, W^) is a 2-dimensional standard Brownian motion. By the strong Markov 
property, K" is a regenerative functional in the sense given in the previous section 
with the stopping times {%"] defined as 



-log(D) - iz,, + ^/Trn^/f'nA"{K"JJTn) 



(26) 



T^' = 0, t" =mfh>T";X^ ^ xl, sup X« > 



with an arbitrarily fixed point (Xp,x") e with x^ < x". Let us take xJJ = x, 
x" = X + a with (x,a) which satisfies Condition 13.21 recall that d„ e C(y,6), so that 
we can find such a pair (x, a) for each n. 

Put s„ = Sg„, n„ - Tig,, and H,, - Hg,, . To verify all the conditions for Theorem l4.6l 
to hold, we use the following more-or-less known identities. The first one is that 

1 



(27) 



n„[g] 



rE 



Jo 



g(X^)dt 
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for all integrable function see e.g., Skorokhod [17). Section 3.1. The second one 
is Kac's moment formula |6l: for given a positive Borel function g, define 



G^„(i/;z) = E« 



■>T(Z) 



r g{X';)G'-\X';;z)dt 



oT(z) 



r i/(xnidi 

Jo 



= fc!Gf^|(y;z) 



recursively for A: e N, where y, z e IR, G^(i/; z) = 1 and 

T(z) = inf{t>0;Xf =z}. 

Then, it holds that 
(28) 

for any y, z e ]R. The third one is that 

Gi{y;z) = 2 r (s„(z) - s„(x))g(x)n„(dx) + 2(s„(z) - s„(y)) T ^(x)n„(dx) 

Jy J-Do 

if y < z, and 

poo 

G^^(y;z) = 2(s„(y)-s„(z)) g(x)n„(dx) + 2 (s„(x) - s„(z))g(x)n„(dx) 

J 1/ Jz 

if y > z. See Skorokhod UtI , Section 3.1 for the details. 

Lemma 5.1. Condition \4.2\ holds. 

Proof: By the strong Markov property and the above identities, 

ml = E14 - tJ] = E«„[t(xJ)] + E:J„[t(x«)] = 2(s„(x») - si(Xo")). 
The result then follows from Condition |3.2| //// 

Lemma 5.2. Condition\4.3\holds. 



Proof: By the Burkholder-Davis-Gundy inequality and the strong Markov 
property, it suffices to show 



SUpEy 



\T{z)f + 



/-T(Z) 

X ' 



\hn{Xl)\dt 



pT(z) 2 

Jo 



< oo 



for (y,z) = (0,Xg), (y,z) = and (y,z) = (x",x[j). We only need to show 



(29) 



SUpEy 



-t(z) 



r (i+(p(x«)^)di 

Jo 



41 



< CO. 



because h„ = T(p^ - E„, 



E„ = Tnn[<p^] = — E" 



pi' 
Jo 



(p(X«)^dt 
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and inf„ m" > by Lemma ISTTl By Condition [XTJ we have 

s'„{w)gk{v)n„{v) < le('^+i)}'.l"l-(4y.+6)K'-'''l) 

for gk{v) = (1 + (piv)^) exp(A:y±|z;|), k e if \v\ > \w\ and vw > 0, where y± = y+ if 
c > and y± = y- otherwise. Hence, by Condition 13.21 there exists a constant C 
(independent of n) such that 

for any m e ]R as long as A: < 3, where z = or z = x". This inequahty implies | |29] l 
with the aid of iglll. //// 



Lemma 5.3. Condition \4.4\ holds. 

Proof: The proof is lengthy so is deferred to Section 6. 
Lemma 5.4. Condition \4.5\ holds. 



nil 



Proof: Note that 



< inf Z„ < sup Z„ < oo 



H6N 



H6N 



by Leimnas l5.1[|5.2l and l5.3l The first two properties are then obvious from l(25l l. To 
see the third, notice that 

T 



V' = 



-E 



J 



(30) 



-E 



r V'„(x«)<p(x;')p(x;' 

Jt" 



)df 



0(e2) 



=l + -fn„[!/.„<pp] + 0(e^) 



and n„[!/;„(j9p] = 0(1), in the light of Lemma I5l2l where 
(31) ^n{y) = 2s'„{y)c„{y) f h„{iu)n„{dw). 

\J — oo 



Here we used the fact that 



r t" li* (x) r ■> 

0= '^dx= !/.„(XndWi+ /z„(Xndt, 
"1 1 1 

which follows from the Ito-Tanaka formula. //// 
Now we are ready to apply Theorem 14.61 In the light of Lemma 3 of Fuka- 
sawa IfTOl and Lemma l5l2l we have A" = 0(e„). Further, by the Ito-Tanaka formula 
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and Lemma [5. 2[ we obtain 



' cp{X^) |p(Xf )dWi + ^1 - p(xpdw2 
= 3E«„ £ cpiXI) |p(X;')dWi + ^1 - p(X«)dWf I (piX'^fdt 



= 3E" 



Jo Jo 



where ipn is defined by | [3T] |. This implies 



Jo 



)p(xr)dwi 



„ 3n„[<ppi/;„] 6a„VT; 
A" = - „ , + 0(e„) = + 0(e„), 



- n„[(p2]3/2r • ■ 

where a - a„is defined as ^ with c - cg^^ - c„/e,„ E - E„ and n = n„. Since 

l^ii OV '0=1 

=c/)(z;l)-a„(/)(z;l)(z2-l) + 0(e2) 

by ( |30] l, we conclude tj"(z) = (^(z; 1)(1 + Ph(z)) + 0{e}j). Hence we obtain a contra- 
diction to 



6. Proof of LemmaII 
Here we prove that the characteristic function V' (u) of 

I, £ K{x^)dt, £ (p{x^i) p(xndwi + ^1 - p(x«)2dw2 

under E"„ satisfies the inequalities of Condition l4.4l By the strong Markov property, 
it suffices to prove the same inequalities for the characteristic function ^"(m) of 

t(x«), /;„(x;')di, I <p(x;') p(x«)dwi + - p(x«)Mwf 1 

under E"„ instead of ^"(m). 

Note that Y" := s„(X") is a local martingale by the Ito-Tanaka formula, so that 
there exists a standard Brownian motion B" such that Y" = B"y») martingale 
representation theorem. Under E"„, Bp = s„(Xg). Note also that 



dY;' = s'„{X1)c„{X'l)dWl = J^{Xl)dWj = -^dwi. 
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where o„(i/) = ^7i„(s„^(i/))/ ^Js'„{s„^{y)). It follows that 

f gix'!)dt = r ^' g{s-HB:))oniB:fdu, 

Jo Jo 

r g(xndwi = r ^(s-i(B«))0„(B:;)dB:; 

Jo Jo 

for every finite stopping time t and locally bounded Borel function g. When 
considering the hitting time t = t(x") of X", we have 

(32) {Y\ = Tn ~ inf{s > 0;B: = s„(xi')}. 

Put y" = s„(x") for z = 0, 1 and y" ^ = s„(2i:g - x"). Notice that by definition, 
inf \y" - y"| > 0, inf \y" - y'\ \ > 0, sup \y" - y\\ < co. 

Lemma 6.1. Let B" be a standard Brownian motion with = y^ and define tn as ( f32l >. 

Let A be a set and g„{-, A) : 'R ^ 'R be a sequence of Borel functions for each A e A with 

(33) sup \gn{v,A)\<°°- 

AeA,«elN,i:ie[y^j,yJ] 

Then there exist positive constants a\ and a2 such that for all A e A and n e N, the 
distribution of 

f"g„{B%A)dt 
Jo 

is infinite divisible with Levy measure L satisfying for all z > 0, 

a? 

L((-co, -z]) V L((z, co)) <ai + —. 

Moreover, 

• if there exists a sequence of intervals I„ c [j/Jj, y'j] such that 

inf|I„|>0, inf g„{v,A)>0, 

«eN /\6A,n6N,Be/„ 

then there exist positive constants a^ and a^ such that 
-a3 + ^<L{{z,co)) 

holds for all z > 0, A e A and n e N, 

• if there exists a sequence of intervals I„ c [yj], y"] such that 

iiif|I„|>0, sup gn{'O,A)<0, 

"^"^ neN,A€A,H€/„ 



then there exist another positive constants ^3 and a^ such that 

_fl3 + ^ <L((oo,z]) 

V|z| 

holds for all z < 0, A e A and n e N. 
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Proof: This can be proved by the same argument as in the proof of Lemma 3 of 
Borisov i3J. //// 

Lemma 6.2. Let {B",B") be a 2-dimensional standard Brownian motion with BJj = y'^ 
and define t„ as (l32] l. Let g„ be a sequence of locally bounded Borel functions with 

sup \g„{v)\ < oo, inf |^„(z;)| > 0. 
Then there exist positive constants a\,a2, fls, fl4 such that the distribution of 

' gniBDdB': 



f 

Jo 



Jo 

is infinite divisible with Levy measure L satisfying 



-fll + — < L((-co, -z]) = L((Z, oo)) < fl3 + — 

z z 



for allz> 0,ne N. 



Proof: Put An - y" - y'^. Let t,7„„ i - 1,. . .,mhe the times at which B" first 
attains the levels y'^ + A„i/m respectively. Put 



Note that /"', z = 1, . . . , m are independent by the strong Markov property. Besides, 
{/,"") !<,<„, is a null array for each n e N since for all e> Q, 

sup P[\C\ >e]< P[Mti/„,N2 > + A,„ 

l<i<m 

which converges to as m ^ oo, where M is a constant, N is a standard normal 
variable independent of B", and 

111 1 
inf B" < v" \ = — - 

-i)/,„<f<T,/,„ *~ -ijj m inf„e]N lyJJ - y'^J ■ 

Hence, /„ is infinite divisible for each n e N. Denoting by L its Levy measure, it 
holds that for every continuity point z > 0, 

m 

(35) limV P[r >z]=L((z,cx,)) 

m—^cc A J 

and for every continuity point z < of L, 

m 

(36) limV P[/r <z] = L((-cx,,z]), 



(34) A,„ = sup P 

l</<m,)ie]N 
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for which see e.g.. Feller fSl, XVII.7. Observe that for z > 0, 

p[r<-z] = p[r>z] 



A„ f A; 

A 

Vmj m 



1 r°° dv 

— Am H I T 

Jmz/{A„ ^fM) ^ + ^ 

where M is a constant. Hence, by L'Hopital's rule and l|35] |, 

L((-oo, -z]) = L((Z, 00)) < fl H 

TIZ 

with a constant a > 0. By the same calculation, we have also 

r // ^^ r // ^^ ^ , i^fwelN A„ VM 

L((-oo, -zj) = L((z, 00)) > -fl + 

7IZ 

with another constant M' > 0. //// 

Lemma 6.3. Let (B", B") fee a 2-dimensional standard Brownian motion with B'^ = y'^ 
and define f„ as (32i . Let Abe a set, gn,i(-. A) be Borel functions for each A e A, g„^2 be 
Borel functions which are absolutely continuous on [y"_y y"] respectively, and g„^3 be Borel 
functions with 

1 

sup |g„,i(^/ A)| V |g„,2(l^)| V \g'„ ^{V)] V \g„,3{v)\ V — < 00. 

A€A,n€K,v€li/"_^,y'l] ' \gn;i\P)\ 

Assume that there exists a sequence of intervals In c [i/[J, y"] with 



inf \ln\ > 

iieN 



such that 



(37) inf gn,i{v,A)>0 or sup gn,i{v,A)<0 

AeA,n6]N,D6/„ AeA,«eN,i.e/„ 

holds. Denote by g„{-; u. A) the characteristic function of } defined as 

J = ui [" g„,i{B^,A)dt + U2 [" g„,2{B\')dB'} + U2 r " g„,3(B'/)dB^, 
Jo Jo Jo 

where u = (mi, U2) e IR^ with \u\ = 1. Then, there exists a constant C e (0, 00) such that 
for every t e R, it holds 



sup \gn{t; u, A)| < Ce 

AeA,jieN,u;|i(|=l 



-VFl/c 
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Proof: Put An = y" - and let nim, i = 1, . . . , m be the times at which B" first 
attains the levels + A„i/m respectively as in the previous proof. Put 



and 7™' = 7™'^ + 7™'^- By the same argument as before, we conclude that } is 
infinitely divisible for each n e N, A e A„ and m e ]R^. We have (|35] | and l l36l l with 
its Levy measure L. Notice that 

(38) gn,2iB'l)dB'l = gnM^V-^ Ci^^D^i 
on the set 

I inf B" > y", 

[T„_i)/„,<f<T,/„, 1 

by the Ito-Tanaka formula. Since, for example, 

p[c >A>p [/,""'' > 2z] - p[r:-' < -z], 

P[C >z]<P [J7'^ > z/l] + PUT'' > z/2], 
there exist positive constants o/, i = 1,2, ... ,6 such that 

(39) - fli + < L((z, oo)) < fl4 + — - + 

yi z Vz z 

for all z > and 

-fll — + — — — < L((-oo,z]) < fl4 + — — + — — — 

VH Z V0 Z 



for all z < by Lemmas l6.1[l6.2l and 

In case we have the first inequality in ll37t , if 

with jSi - mf,\eA,n€N,veI„ gn.iiv. A) and jS2 = SUp„gj^ pgj^^ l^,',,2(^)l' t^*^" 

inf Mig„,i(u; A) - M2^;,2(u)/2 

> Vl - I"2P/51 - |M2liS2/2 > iSl/4 > 0, 

so that by Lemma |6A] and i 



limV P[7r'>z]>-fli + ^ 
for all z > 0, where fl,, i = 1,2 are positive constants. In addition, we have 

Pir >A> PUT > > z] > ^PUT'' > z\ 
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for all z > 0. Hence, when \u2\ < [io, there exist another constants a'., i = 1,2 such 
that 

(40) - + 4 ^ CO)) < + ^ + ^ 

for all z > 0. Now, note that by the Levy-Khinchin expression, there exists a 
constant ct^ > such that 

Re log(^„(i; u, A)) = -if sin^{zt/2)L{dz). 

Jr 

Take Zi > such that z i-^ sin^(z) is increasing on [0,zi/2]. Then, observe that for 
sufficiently small zq e (0, zi), it holds that 

ay y/zo - as/ ^ fl2/ + fls/ ^ ^ 
fl6/zi fls/zo-fle/zi 
Fix such a point Zp and take jSs such that 

ay ^/z^- as/ ^/z^ ail ^[2^ + a^l ^fzl 

Then we have for the case that IM2I Vf > jSa, 
sin2(zt/2)L(dz) 



I 

Jr 



> I sin2(zt/2)L(dz) 
(41) >sin2(zo/2)L((zo/|fUi/|f|]) 



> sin2(zo/2) ((^ - ^) \uM - + ^) ^ - 

> ^^/c - log(C) 

for sufficiently large constant C by | |39] l. Further, by (|40l l, we have for the case that 
IW2I Vf < jSs, 

sin2(zo/2)L((zo/|4zi/|i|]) 

> sin2(zo/2) |(^ -^]^\- ^\u2\\t\ - a\ - aX 

> ^^/c - log(C). 

The same conclusion is obtained also in the case that we have the second in- 
equality instead of the first in | |37t . For example, we define /3i alternatively as 

h = - supAeA,„eN,i,e/„ gn,i{v,A) and observe 

m 

lim Y < z] > -5i + ''^ 
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for all 2 < with positive constants ai, 02 when \u2\ < jSo- Then use 
PUT < z] > PUT < z; <z]> ^PUT'' < z] 



to obtain 



-a[ + —p= < L((co,2]) < fl4 + — p= + 



for all 2 < 0. The rest is a straightforward translation. //// 
Now we are ready to prove Lemma 15.21 By Petrov's lemma (see Petrov |16|, 
p. 10), it suffices to prove that there exists a constant C e (0, 00) such that 

|^"(m)| < Ce-I"l/^, 

where is what was defined at the beginning of this section. Put 

gn,iiv,A) = (Ai + A2h„{s~^{v)))a„(v), 
gn,2{v) = (p{s-\v))p{s~\v))a„{v), 

and 

M(A) = sup |g„,i(z;, A))| V \g„,2{v)\ V 2(1^)! V |^„,3(^')l V — ^ 

n6N,i>6[i/;^j,y;] \gn,3\,<J)\ 

for A = (Ai, A2) e S, the 1-dimensional unit sphere. It is not difficult to see that 
for all A e S, we have M(A) < 00 and that there exists a sequence of intervals 
J„(A)c[y^',y«]with 



such that 



inf \In{A)\ > 

neN 



inf gn\{^,A)>0 or sup gn\{^,A)<0 



holds. If the first inequality holds for A - Aq, put 



and 



m(Ao) = inf gnAv, Ao) 

n€N,v€l„{Ao) 



A(Ao) = A e S; inf g„,i{v, A) > m(Ao)/2,M(A) < 2M(Ao) } . 

n€K,v€l„{Ats) 



If the second inequality holds for A = Aq, put 



and 



m{Ao)= sup ^„,i(i7,Ao) 

I16]N,P€/„(Ao) 



A(Ao) = A e S; sup g„,i{v, A) < m(Ao)/2,M(A) < 2M(Ao) 

I n6lN,ti6l„(Ao) 
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Now, notice that A(Ao), Ao e S is an open covering of S, so that it has a finite 
subcovering A(Ai), ■ ■ ■ , A(Aj). For each A^, we can apply Lemma [63l with A = A(Ay) 
to obtain that there exists Cy > such that 

sup \W"ituiAi, tuiAi, tu2)\ < Cye-I'l^'^' 

AeA(Ay),w€S,ne]N 

for all t e ]R. Since / < oo, we conclude that there exists C > such that 
sup |4'"(fMiAi, h/iAz, tu2)\ < Ce-I'l/^ 

AeS,ueS,«6N 

for all t e ]R, which completes the proof of Lemma [5.2| 
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